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ABSTRACT 

We present a global study of low mass, young stellar object (YSO) surface densities 
(S) in nearby (< 500 pc) star forming regions based on a comprehensive collection of 
Spitzer Space Telescope surveys. We show that the distribution of YSO surface densi- 
ties in the solar neighbourhood is a smooth distribution, being adequately described 
by a lognormal function from a few to 10 3 YSOs per pc 2 , with a peak at ~ 22 stars 
pc~ 2 and a dispersion of o"i O g 10 s ~ 0.85. We do not find evidence for multiple discrete 
modes of star-formation (e.g. clustered and distributed). Comparing the observed sur- 
face density distribution to previously reported surface density threshold definitions 
0^ ' of clusters, we find that the fraction of stars in clusters is crucially dependent on the 

adopted definitions, ranging from 40 to 90%. However, we find that only a low fraction 
(< 26%) of stars are formed in dense environments where their formation/evolution 
(along with their circumstellar disks and/or planets) may be affected by the close 
proximity of their low-mass neighbours. 
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1 INTRODUCTION 

It is often stated that most if not all stars form in stellar 
clusters. This view is based largely on near-infrared (NIR) 
studies of star-form ing (SF) regions within several hundre d 
parsecs of the Sun (|Lada fc Ladall2003l ; iPorras et ai1l2003l) . 
However, adding high-resolution mid-infrared (MIR) data 
to the NIR makes YSO identification more robust and less 
likely to be contaminated by field stars, which leads to bet- 
ter tracing of YSO surface densities. This means that with 
the NIR alone, there were large uncertainties in the num- 
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ber of stars at lo w values of YSO surface densities (Syso) 
l|Carpenterll2000h . 

The spatial distribution of forming stars, i.e. do they 
form in clusters, is important for two main reasons. The 
first is that dense environments can affect the evolution of 
the young stars as well as alter th eir disk and planet forma- 
tion/evolution (All en et al.l 120071 ). The second is to locate 
the progenitor population of open clusters and to determine 
why such a low fraction of the Galactic stellar population 
is observed in clusters. Are there multiple discrete modes, 
such as clustered and distributed, in the star-formation pro- 
cess that manifest themselves as peaks in a surface den- 
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Name 


YSO Number 


Distance 
pc 


Reference 


Auriga 


138(172) 


300 


1 


Ccphcus I 


34(46) 


280 


1 


Ccphcus III 


44(52) 


280 


1 


Cepheus V 


19(19) 


280 


1 


Chameleon I 


67(93) 


200 


1 


Corona Australia 


27(45) 


130 


1 


Lupus III 


43(79) 


150 


2 


Ophiuchus 


199(297) 


125 


2 


Orion* 


2696(3352) 


111 


3 


Perseus 


280(387) 


250 


2 


Serpens 


179(262) 


415 


2 


Taurus 


131(249) 


137 


4 



Table 1. The Spitzer surveys used in the present work includes 12 
star forming regions with 3857 YSOs. The numbers in brackets 
refer to the total number of sources in the catalogues for each 
region, while the number before the brackets is the number used 
in the present analysis. The difference is due to the application of 
the absolute magnitude cuts as well as the elimination of class III 
YSOs from the sample. The sources for these SF regions are the 1) 
GB survey, 2) c2d survey, 3) Orion survey and 4) Taurus survey. 



"ONC is excluded, see §[2] 



. .Strom et all 1 19931 : ICarpenterl [20001: 
Wang et al.ll2009h ? 

of the Spitzer Space Telescope 



sity distribution (e .g 
IWeidner et ai1l2004l :h 

With the lau nch 
I Werner et alj|2004l ) we are now able to differentiate YSOs 
and contaminating sources based on colour information and 
hence can study the distribution of YSOs independently 
of the surface densities. La rge field-of-view (FoV) Spitzer 
obser vations of SF regions l|Allen et all 120071 : H vans et alj 
l2009h found that YSOs extend well beyond the densest 
groups in their environment and continue throughout. We 
combine several Spitzer surveys that cover nearly all the SF 
regions within 500 pc of the Sun. A list of the regions and 
their properties is given in Table [T] Note that with only the 
local SF environments being considered, we are not sampling 
massive star forming regions that are found beyond 500 pc. 

Using the comprehensive collection of Eyso we inves- 
tigate what fraction of YSOs are found in dense clusters. 
We define dense clusters as regions where YSOs are af- 
fected by their neighbours in sufficiently short timescales 
of < 10 5 yr, such that its surface densi ties exceed ~ 
200 YSO pc -2 (see iGutermuth et al l [20051 ). We also re- 
view what surface densities are required to identify "clus - 
ters" accordin g to de fi nitions provi d ed by [Carpenten fe OQO): 
lLada fc Ladal J2003I); lAllen et all (|2007l ); Ijprgensen et all 



( 2008); IGutermuth et al.1 (|2009l ) in §[5] In this paper we will 



investigate 1) whether there is evidence for multi-modality 
in the surface densities of YSOs, 2) what fraction of stars 
form in dense clusters in the local neighbourhood and 3) 
how relevant the various cluster definitions are. 



2 OBSERVATIONS & DATA 

Multiple Spitzer surveys were used to generate a compre- 
hensive and statistically significant dataset to investigate 
the spatial surface density properties of forming stars in the 
solar neighbourhood. The surveys are the Gould's Belt (GB) 
survey (Allen et al. in prep.), Orion survey (Megeath et al. 



in prep.), Cores to Disk s (c2d) survey (lEvans et al.ll2003l ). 
and the Taurus survey (|Rebull et aUl2010l '). The GB and 
Orion catalogs have not been publicly released yet. We have 
more than 7000 YSO detections in the combined catalogs at 
distances between 100 to 500 pc. 

Spitzer data are necessary for this study as low Eyso 
can be differentiated from field star populations, unlike NIR 
observations where field star contamination can be problem- 
atic. The YSO population that we have collected represents 
a global view of the low-mass star-forming region in the lo- 
cal neighbourhood from low to high surface densities. These 
Spitzer surveys combined represent the most complete cen- 
sus of star formation within 500 pc of the Sun available to 
date. 

In order to homogenise the data from the surveys we ac- 
counted for distance effects on photometry, namely we limit 
the absolute magnitude range used for individual sources 
to that of the faintest YSO detectable in the furthest SF 
region and the brightest in the nearest SF region. The ab- 
solute magnitude limit used for the 500 pc data collection 
is ^ M3.6 U m ^ 5.91, based on Orion at a dis tance of 
414 pc i|Menten et al.ll2007l : iMavne fc Navlorll2008h for the 
faint sources and Ophiuchus at 125 pc for the bright sources. 
This reduces the number of YSOs we can use, but it miti- 
gates detection biases introduced for SF regions at different 
distances. 

The GB and c2d surveys classify Y SOs using spectral 
indices (|Ladal Il987l : iGreene et al]|l994l ). The Taurus and 
Orion YSOs are c l assified by using colour-colour diagrams 
IIAllen et al. | 120041 : iMegeath et al. M2004I : IGutermuth et al.l 
120051 . 120091 ) What fraction of the YSOs are diskless, gen- 
erally classified as Class III, and h e nce n ot identifiable in 
the IR? Based on iHernandez et all (|2007l ) we assume that 
65% of the YSOs have disks. We corrected the stellar surface 
densities of the data for the missing fraction of 35%. 

Orion, which offers the largest range of stellar surface 
densities and hosts the most massive stars of the SF regions 
considered in this paper, had to be treated separately from 
the other surveys. The ONC, in particular the Trapezium re- 
gion, has two Spitzer based issues: Stellar surface densities 
that exceed Spitzer' s spatial resolution and the extremely 
bright nebulosity that diminishes effective sensitivity con- 
siderably. The bright nebulosity introduces errors for YSO 
identification since the PAH emission outshines lower mass 
YSOs and introduces large errors in the photometry. To 
compensate for the complex incompleteness we removed all 
YSOs centred on Ol Orionis within a radius of 0.56 pc (4.7'). 
To correct for missing YSOs from the removed region, we es- 
timated t hat the mass remove d was ~ 25% of the total Orion 
complex l|Getman et al.|[2005l ). Excluding the ONC from our 
analysis does not significantly change the presented cumula- 
tive distribution of surface densities presented in this paper. 
If we were able to observe all the members in the ONC based 
on the ~ 25% of mass we estimated to be missing, the av- 
erage LonC < 1000 YSOs/pc 2 . This surface density regime 
goes beyond the scope of values we are presently consider- 
ing. Hence we are not sensitive to the extreme high E tail 
end of the ONC distribution. 

Spitzer is not completely free of co ntamination when 
identifying YSOs, i.e . AGBs/Be s tars ([Cieza et al] 
Robitaille et alj|2008f) and galaxies (IGutermuth et ail '. 



2010; 



2008; 



Evans et al.ll2009l ). lOliveira et all (|2009l ) found that ~ 25% 
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Figure 1. Top panel: The surface density distribution of the to- 
tal sample of YSOs in the solar neighbourhood used in this work 
(black). A lognormal function with a peak at ~ 22 YSOs/pc 2 and 
a dispersion <ri ogloS = 0.85 is shown as a dashed (red) line. Bot- 
tom panel: The same as the top panel but now broken into the 
three respective surveys. Note that Orion dominates the number 
statistics. 



of the identified YSOs in the c2d Serpens catalog are AGBs, 
which is likely an isolated worst case scenario as Serpens is 
the field closest to the Galactic plane in our compilation of 
SF re gions. Two of the twenty contaminants lOliveira et al.l 
(2009) identified are Class lis and the rest of the contam- 
inants are Class Ills. We only consider Class I /II objects, 
where the AGB contamination is < 10%, and remove all 
Class Ills. The flat spectrum sources are grouped with Class 
I objects. Between the methods used to identify YSOs in 
the c2d, GB, Tau rus and Orion data, which are the c2d 
jEvans et al.ll2009T ) and Gutermuth et al (2008, 2009) meth - 
ods, the selection discrepancy is < 5% l|Rebullet alJl20ld ). 
By selection discrepancy we mean the agreement that an 
object is or is not a YSO (Class I/II). 

Extra-galactic background contamination for YSO MIR 
identification is well studied. For the c2d and GB catalogs, 
which use the same data-reduction pipeline, I Evans et al.l 
(2009) found that background galaxies contaminate ^ 
5% of the YSOs. Si milarly, YSOs identified via the 
iGutermuth et all (|2009r i method for Orion is < 1%. For Tau- 
rus t he expected contamination rate is ^ 5% |Rebull et al.l 
l2010h . 



3 Eyso DISTRIBUTIONS 

Our primary tool for analysing the surface densities is com- 
puting the local observed surface density of YSOs centred on 
each YSO's position, where Eyso = (N — 1)/(ttD%) and N is 
the Nth nearest neigh bour, and Djy is the proj ected distance 
to that neighbour (see lCasertano fc HutHl98a ). Throughout 
this work we will adopt N = 7, although we note that all 
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Figure 2. (a) The cumulative fraction of surface densities for the 
GB+Taurus, c2d, and Orion surveys. Each SF region included in 
the distributions has N(YSOs) ^ 10 and a sufficient field-of-view 
to properly calculate stellar surface densities. The Orion survey 
stops at 73% for the cumulative fraction since the ONC is ex- 
cluded. We adopt a 65% disk fraction for all of the SF regions. 
We normalised each curve by the number of YSOs in each survey, 
(b) With the GB+Taurus, c2d, and Orion surveys combined we 
see Class I & II distributions having similar profiles with a small 
offset in density, showing that we are likely seeing the primor- 
dial distribution of the YSOs. (c) With all of the Spitzer surveys 
combined we compare several cluster definitions. The vertical grey 
lines from left to right are Lada & Lada (2003), Megeath et al. 
(in prep.), j0rgensen et al. (2008), Carpenter (2000), and Guter- 
muth et al. (2009) stellar density requirements for clusters. These 
values correspond to 3, 10, 20, 32, and 60 YSOs pc -2 and inter- 
sect the corrected cumulative distribution profile, implying that 
87%, 73%, 62%, 55%, and 43% of stars form in clusters, respec- 
tively. The percentages correlate to what fraction of stars form 
in "clusters" based on the various definitions. The black vertical 
line is for a dense cluster where E ^ 200 YSOs/pc 2 . The fraction 
of YSOs in a dense cluster is < 26%. 



results have been tested for N = 4 — 22 and no significant 
differences were found. 

Figure [JJ shows the surface density distribution of all 
YSOs in our sample, corrected for the diskless fraction. Ad- 
ditionally, we show a lognormal fit to the data as a dashed 
red line (see § [4} . The over-prediction of the lognormal at 
high Eyso compared to the observations is most likely due 
to the exclusion of the ONC and surrounding area (see §[2]). 
The bottom panel of Fig. [JJ shows the surface density distri- 
bution for each of the three surveys separately. 

In order to see the fraction of YSOs above a given E 
threshold, we show the combined Eyso distribution (shown 
as a cumulative fraction normalised to the number in each 
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combined survey) for the three surveys used in this study 
in Fig. [2^,. Note that the GB/Taurus distribution lies to the 
left of the c2d survey. This is simply due to the GB/Taurus 
focussing on lower density regions than c2d. The cumula- 
tive distribution for the Orion survey only reaches 0.73 in 
Fig. [2^ and 0.81 in Fig. [2};, where all the surveys have been 
combined, due to the exclusion of the ONC. In Fig. [2J; we 
show the cumulative distribution of all YSOs included in our 
survey, while in Fig. [Sp we split the survey into class I and 
class II objects. 



4 RESULTS 

It has been long assumed that two distinct modes of star 
formation exist for YSOs, 'clustered' and 'distributed' (e.g. 
iGomez et al.ll 19931 ; ICarpenterll2000l ; lLada fc Ladall2003h . but 
the notion has been questio ned after Spitzer results hinted 
otherwise (I Allen et al.ll2007h . If there are indeed two modes, 
then we would expect to see a bi- or multi-modal profile in 
cumulative surface density distribution plots such as Figs. 
[1] [2^i, & [2};. Instead we see smooth and featureless dis- 
tributions from the low to high stellar surface densities for 
the c2d, GB, Taurus, and Orion surveys. We find that the 
Eyso distribution of low-mass stars in the solar neighbour- 
hood can be well described by a lognormal function, as seen 
in Fig. [1] with a peak at ~ 22 YSOs/pc 2 and a dispersion 
ciog 10 s = 0.85. 

The spatial distribution of the YSOs in these SF regions 
is expected to be close to primordial since their YSOs, in par- 
ticular Class Is and Class lis , are ^ 2 Myr old l|Haisch et al.l 
l200ll ; lHernandez et al.ll2007f ). In order to place stricter con- 
straints on this, we now split the complete sample into Class 
I and II objects, which can be roughly attributed to an age 
sequence. The cumulative E distributions of Class I and II 
YSOs are shown in Fig. [5p. We see that the two distribu- 
tions have similar smooth density spectra, however they are 
slightly offset. The E of the Class I/II objects are calculated 
by finding a YSO's Nth nearest YSO. Once this is done 
for the YSOs we separate the Class I/II objects. E is cal- 
culated this way since Class Is and Class lis are not alway s 
spatially distinct from one another dGutermuth et al] [2009). 
Class lis are known to be slightly more dispersed than Class 
Is in high density regions l|Gutermuth et al]|2009l ) reflecting 
early dynamical evolution. However, the similar distribution 
between these classes leads us to conclude that the distribu- 
tion of observed E is mainly primordial in nature. 



5 CLUSTER IDENTIFICATION 

The definitions of what defines a cluster vary widely as we 
have limited knowledge about YSO membership other than 
their projected two dimensional spatial distri butions. Some 
defin itions have a physical motivation (e.g. lLada fc Ladal 
2003) while the others are gen erally empirically -derived 
from the data being considered (| Allen et al.1 120071 ). When 
applied to a uniform dataset like ours, differing choices of 
a surface density threshold returned different "clustered 
fractions", as summarised below. 



I Carpenter! (|2000h : Clusters in SF regions are identi- 
fied by using stellar density maps in the K s band. The 
density maps are field-star background-subtracted (galactic 
coordinate dependent) based on semi-empirical models. 
Clusters are identified as 2a over-densities and defined as 
regions with 6a over-densities (with total number of mem- 
bers taken as the number of sources above the 4a threshold) 
with respect to the local background. Carpenter's cluster 
Eyso ranged from 20 to 67 YSOs pc -2 with a median of 32 
YSOs pc -2 . Considering the median, 55% of the YSOs are 
contained in clusters. 

lLada fc Ladal ((2003): A physically related group of 
stars, called an embedded cluster, that is 1) partially or 
fully enshrouded in interstellar gas and dust, 2) has ^ 35 
YSOs and 3) a stellar-mass volume density of 1.0 Mq 
pc -3 or greater such that its evaporation time exceeds 10 s 
years. In surface density, rather than volume density, the 
number of YSOs pc~ 2 necessary for "cluster" is ~ 3 (see 
|j0rgensen et al. 1 120081 ). The authors estimated that 80-90% 
of the YSOs are in embedded clusters, which is found to be 
in agreement with our Spitzer data. 

|j0rgensen et all (|2008h : Building upon the lLada fc Ladal 
(2003) definition of an embedded cluster, j0rgensen et al. 
define a cluster as being "loose", which is the same as an 
embedded cluster, and a "tight" cluster. A tight cluster 
requires a stellar-mass volume density of ^ 25 Mq pc -3 
and > 35 YSOs, which implies that 62% of the YSOs from 
our data are contai ned in such c l usters . This finding is close 
to 54% as found in lEvans et al] ([20091 ). 

iGutermuth et al.l (|200 9): This method employs the 
minimal spanning tree (MST) algorithm to define cluster 
cores by isolating the densest parts of larger scale over- 
densities. The MST is a network of lines that connects a 
set of points, has no closed loops, and the set of edges add 
up to the shortest total length possible between all points. 
After determining a cutoff length for the MST collection, 
YSOs can be separated into two populations: clustered 
and distributed. The authors found that the clusters from 
this analysis range between 0.64 and 78 YSOs pc -2 with 
a median of 60 YSOs pc -2 . Roughly 43% of the YSOs are 
found in a median core clusters. 

Megeath et al. in prep.: A cluster is a set of con- 
tiguous objects which have nearest neighbour densities 
10 YSO pc" 2 . The 10 YSOs p c ~ 2 is similar to the cluster 
definition given in lAllen et all <|2007l ) and motivated by a 
comparison of the Orion (Megeath et al. in prep.) and the 
Taurus molecular clouds. The Taurus and other similar 
dark clouds, i.e. Chameleon and Lupus, have most of their 
objects at densities below 10 pc -2 , while Orion and other 
clouds with clusters have 70-80% above this threshold. 
Applying the 10 YSOs pc -2 definition to our dataset results 
in 73% of YSOs being in clusters. 

In Fig. [2}; we show five vertical grey lines that refer 
to the defined densities required for a collection of YSOs 
to be considered "clustered" (Lada & Lada 2003; Megeath 
et al. in prep.; j0rgensen et al. 2008; Carpenter 2000, 
Gutermuth et al. 2009). The vertical lines fall on the same 
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featureless slope and do not correspond to any preferred 
density. The black vert ical line, which correspo nds to dense 
clusters (as defined in iGutermuth et al.l 12009 ). shows that 
< 26% of YSOs are formed in environments where they 
(along with their disks and planets) are likely to interact 
with their neighbours. 



6 DISCUSSION AND CONCLUSIONS 

We have compared our global surface density distribution 
with previously reported definitions of clusters (discussed 
m § 0), an d nn d that the fraction of stars in the solar 
neighbourhood forming in clusters is crucially dependent 
on the adopt e d defi nitions (ranging from ~ 40 to 90%). 
lLada fc Ladal (|2003h used a physically motivated definition 
of clusters, and their adopted low surface density of ~ 3 
YSO pc -2 encompasses nearly all star formation in the so- 
lar neighbourhood. However, only a small fraction (< 26%) 
of stars form in dense clusters where their formation and/or 
evolution is expected to be influenced by their surroundings. 

We conclude that stars form in a broad and smooth 
spectrum of surface densities and do not find evidence for 
discrete modes of star formation in the £ of low mass YSOs 
forming in the solar neighbourhood. Only a small fraction 
of YSOs form in dense clusters where nearby YSO members 
affect its disk/planets evolution. The observed lognormal 
surface density distribution is consistent with predictions 
of hierarchically structured star-formation, wher e the struc- 
ture comes from the MC hierarchical structure (Elmegreen 
120021 . 1 2008). By hierarchical structure we mean a smoothly 
varying non-uniform distribution of densities, where denser 
subareas are nested within larger, less dense areas l|Scald 
ll985l ; lElmegreen et al.ll2006l ; lBastian et alj2007l ). Star form- 
ing environments provide the initial conditions from which 
star clusters may eventually form, albeit rarely. Since the 
probability density function of molecular gas varies with en- 
vironment, as does the tidal field experienced by the SF re- 
gion, it is likely that the fraction of YS Os ending up in b ound 
star clusters varies with environment (|Elmegreenll2008T ) and 
the observed Eyso is not universal. Hence, in a future study 
we will extend this work out to 2 kpc, which includes high- 
mass star-forming regions and more extreme environments 
that may show different results than what we see for the 
solar neighbourhood. 
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